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In this paper we study numerically and experimentally the possibility of using metallic photonic 
crystals (PCs) of different geometries (log-piles, direct and inverse opals) as visible light sources. 
It is found that by tuning geometrical parameters of a direct opal PC one can achieve substantial 
reduction of the emissivity in the infrared along with its increase in the visible. We take into account 
disorder of the PC elements in their sizes and positions, and get quantitative agreement between 
the numerical and experimental results. We analyze the influence of known temperature-resistant 
refractory host materials necessary for fixing the PC elements, and find that PC effects become 
completely destroyed at high temperatures due to the host absorption. Therefore, creating PC- 
based visible light sources requires that low-absorbing refractory materials for embedding medium 
be found. 
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1. INTRODUCTION 

Photonic crystals (PCs) are artificial structures char- 
acterized by periodical variation of the dielectric func- 
tion in space. Their optical properties differ from the 
optical properties of continuous media: electromagnetic 
waves in the PC have a band structure similar to that 
of electrons in a solid crystal I n particular, the abil- 
ity of a PC to modify the thermal emission has drawn 
much attention recently [2T-fl7l]. Generally, the thermal 
emission from a medium at a certain wavelength is de- 
termined by an interplay between the photonic density of 
states and the energy transport velocity 0, Q . Thus the 
thermal emission from a PC exhibits strong wavelength 
dependence, for example, it may be suppressed for wave- 
lengths corresponding to band gaps, where the photonic 
density of states is zero, or stop bands where it is zero for 
certain wave directions. By changing geometric param- 
eters of a PC, one can tune the high or low emissivity 
wavelength ranges to the desired parts of the spectrum. 
This effect has been studied theoretically and experimen- 
tally for different PC geometries, such as log-piles 0-[l|, 
opals [13, |Tj| , inverse opals flgl . pyra mids (l3j , inverted 
pyramids 14], arrays of holes [15l - Fl7T | . etc. 

One of the most important applications for PC- 
modified thermal emission is thermophotovoltaics Q. 
There the spectrum of an emitter should be adjusted 
to the photovoltaic cell sensitive wavelength range in 
the infrared (IR) [l8[. Additionally to PC bandstruc- 
ture effects, possible further ways to enhance efficiency 
of thermophotovoltaic systems can be based on plas- 
monic near-field emitters [ill [2(| , negative index meta- 
materials [U [22J, intersubband transitions in multiple 
quantum wells tailored with photonic-crystal resonant 
effects 23]. Photovoltaic cell by itself can be designed 
as a PC with improved solar power conversion efficiency 
which is achieved due to enhanced absorption [24j and 
applying curved pn-j unction [25j . 



Recently it was proposed to use PCs as a source of 
visible light 4, l2a - [2q| . Suppression of the IR emission of 
such sources may increase their efficiency in comparison 
with conventional incandescent lamps which emit mostly 
(about 95%) in the IR. Some candidate PC geometries 
have been studied theoretically so far. It was shown that 
it is possible to achieve high emissivity in the visible 
range using one-dimensional tungsten photonic crystal, 
however, emission in the near IR appears to be high as 
well [27[. Emission in the IR could be suppressed using 
iridium inverted square-spiral photonic crystal [28[, but 
fabrication and use of such a structure as a light source 
is impractical due to high costs of iridium. 

In this work we focus on practically achievable PC- 
based visible light emitters, namely the three- and two- 
dimensional PC structures with tungsten emitting ele- 
ments. Tungsten is widely used as a traditional light 
emitting material for incandescent lamps due to its high 
melting temperature (3695-fT) and high selectivity of 
emission in the visible range. We analyze a number of 
possible three-dimensional tungsten PC geometries (log- 
piles, direct and inverse opals), perform numerical simu- 
lations, and support the obtained results experimentally. 

The rest of the paper is organized as follows. In sec- 
tion II, we describe the method for the emissivity cal- 
culations and briefly discuss its accuracy. In section III, 
we show that log-piles have a poor emission selectivity in 
the visible and cannot be used for designing efficient light 
sources. In section IV, we find that by choosing appropri- 
ate geometrical parameters of a direct opal PC, one can 
achieve low emissivity in the IR and high emissivity in 
the visible. In section V, we study the effect of position 
and size disorder of the PC elements and present com- 
parison with experimental results. In section VI, we take 
optical properties of the refractory host medium neces- 
sary for fixing the direct opal PC elements into account. 
In the final section VII, we summarize our results. 
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2. PC EMISSIVITY CALCULATION METHOD 

According to Kirchhoff 's law of thermal radiation, the 
emissivity of a structure is equal to its absorption coef- 
ficient in thermodynamic equilibrium Q. This law may 
be used for obtaining the emissivity in a numerical ex- 
periment by calculating absorption instead of simulat- 
ing radiation from a structure. Direct numerical simula- 
tion of a radiating body is possible by various techniques, 
for example, by using stochastic Langevin electrodynam- 
ics [29|, [3(| or by truncating the natural modal expan- 
sion of electromagnetic fields inside a weakly absorbing 
PC [3l|. In both cases rigorous checks of Kirchhoff 's law 
were performed to prove the consistency of the proposed 
models for radiation in thermodynamic equilibrium. 

In the present work we apply the finite-difference time- 
domain (FDTD) method 32] to obtain the absorption co- 
efficient of a structure numerically as a function of wave- 
length. The simulation software was especially created 
for this problem based on our parallel Electromagnetic 
Template Library [33[ • The use of FDTD has a number 
of advantages for simulating PC-like structures. It is able 
to handle purely periodic structures with 2D or 3D pe- 
riodicity as well as quasi-periodic geometries with large 
number of randomly disordered unit cells. For large sim- 
ulations the method is parallelized via domain decompo- 
sition and scales almost linearly with the number of used 
CPUs 33]. Absorbing materials are simulated by fitting 
the experimental complex dielectric function by a num- 
ber of Drude-Lorentz [3 21 ] or modified Lorentz (34| terms. 
There is no restriction on the form of scatterers in FDTD, 
however a special treatment of fine structures such as thin 
material interfaces is sometimes required (35j . 

For numerical experiment we use a conventional FDTD 
calculation scheme in which a simulated PC structure, 
defined by the complex dielectric function of the ma- 
terial is placed in the simulation volume. We use the 
orthogonal uniform Yee mesh for the whole volume. A 
plane wave pulse is generated on a surface outside the 
structure with the help of the Total Field/Scattered Field 
technique and propagates through it (Fig.[TJ inset). Gen- 
erally, a PC consisting of several layers stacked in the Z- 
direction is placed inside a computational volume. The 
plane wave pulse impinges on the PC at given direction 
characterized by the incident angles 9 and <f). We sim- 
ulate oblique incidence of a plane wave on the PC by 
the iterative FDTD algorithm [36[. The computational 
volume is terminated from both sides in the Z direction 
by absorbing boundary conditions - perfectly matched 
layers (PMLs) [32]. They simulate withdrawal of the 
transmitted and reflected waves to the infinity. Peri- 
odic boundary conditions are applied along the X- and 
Y-borders of the computational volume. A residual nu- 
merical reflection from the PMLs is reduced by applying 
additional back absorbing layers technique [33]. Dur- 
ing the numerical experiment, the reflected and trans- 
mitted waves are recorded behind and in front of the 
structure, transformed to the frequency domain and nor- 



malized to the incident spectrum. Ultimately, the trans- 
mission 7{ui,6,4>), reflection 31(lu , 9 , (f>) , and absorption 
A(u},0,(p) = 1 — 7(uj 7 9,(f>) — "5L{u),6,4>) coefficients are 
calculated. The emissivity is then obtained by applying 
Kirchhoff's law as £(w, 9, <fi) — A(lu, 9, <fi). 

Throughout this paper the following FDTD parame- 
ters are used: a) for the case of log-pile PCs, the mesh 
resolution is 100 steps per the PC unit cell (Ar = 42nm 
for the log-pile with the lattice constant a = 4.2/xm and 
Ar = 5nm - for a = 500nm, see Section III), time 
step cAt = 0.5Ar, the numerical experiment duration is 
Nt = 40000 iterations; b) for the case of opal PCs (typ- 
ical PC parameters: lattice constant a ~ 500 nm, sphere 
radius R > 50nm), Ar = 2.5nm, time step cAt — 0.5Ar, 
the numerical experiment duration is Nt — 20000 itera- 
tions. 

FDTD is a well established method for simulating wave 
propagation. The accuracy of any particular simulation 
should however be checked by the analysis of its con- 
vergence with respect to the decreasing mesh step. In 
cases of poor convergence, special FDTD techniques, for 
example subpixel smoothing for interfaces or additional 
PML layers may be required. Convergence and accu- 
racy of FDTD calculations based on EMTL have been 
demonstrated in our previous papers [H, • In par- 

ticular, we considered a general case of an obliquely inci- 
dent plane wave on a metallic opal PC [36j . and demon- 
strated an excellent agreement of the FDTD results with 
the results obtained by the multiple scattering formal- 
ism [H, [39j]. Also, the opal PC case requires a special 
care as the spherical shape of the PC elements should be 
adequately represented on the rectangular FDTD grid. 
In our previous work we presented a subpixel smooth- 
ing method based on the inverse dielectric tensor ap- 
proach [3{| . This method has been shown to significantly 
improve the accuracy of the FDTD simulations for arbi- 
trary shaped metal scatterers and is used in the present 
work for the opal PCs simulations. 

In the present work tungsten is considered as the emit- 
ting material for PC light sources. The data on the tung- 
sten dielectric function e(w) in Drude-Lorentz form is 
taken from [401 ]. Our numerical results are obtained for 
the temperature T = 2400-ftf (if not stated otherwise) 
which is the typical working temperature of a tungsten 
filament in incandescent lamps. 

Energy flow from a PC layer with the emissivity 
£(w,#, (j), T) operating at an ambient temperature T 
within a given frequency range u>± < to < L02 is calcu- 
lated as 

J(T) = J dto J d<j> J Z(uj,9,(t>,T)u(uj 1 T)cos9sm9d9, 

ui 

(1) 

where u(ui, T) is the blackbody spectral radiation density 

hio 3 1 
T) = 4^c 2 exp^/fcT)-!' (2) 
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The efficiency of light emission E is obtained as the ratio 
of the emitted energy flow in the visible J v ; s to the total 
energy flow J tot . 

Note that typical thickness of a tungsten glower (~ 
0.2mm) is much larger than emission wavelength (both 
in visible and IR) . Therefore to calculate the efficiency of 
the glower one can simulate it as a semi-infinite tung- 
sten layer. For such a layer, we obtain an efficiency 
of £\v = 5% and emitted energy flow in the visible 
J vis W = 30kWt/m 2 (at T = 2400X). We use this val- 
ues as references characterizing a standard incandescent 
lamp. Our purpose is to find such a tungsten PC geom- 
etry which yields the efficiency E several times greater 
than Ey/, keeping J V1S of the same order as J v i s ,w at the 
same time, so as to provide sufficient emitted energy in 
the visible. 
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Figure 1: Normal incidence reflection and absorption spectra 
of a log-pile with number of layers N = 4, rod spacing a = 
4.2/1171, rod width w = 1.2/im, and rod height h = 2.4/im. 
High reflectance region at A > 8pm correspond to the IR band 
gap. Sharp peak appears at the IR band gap edge. Inset: 
scheme of the FDTD simulation geometry; 1 - generating 
(Total Field / Scattered Field) border; 2, 2' - detector arrays 
for reflected and transmitted signals. 



Metallic log-pile structures have drawn a special atten- 
tion due to their possible use in thermophotovoltaics [|| 
as IR emitters. Particularly, tungsten log-piles provid- 
ing high emission selectivity in the near IR range were 
demonstrated 5-8]. These structures are characterized 
by a photonic band gap in the IR range (we refer to it 
as the IR band gap later on), and a stop band in TX 
direction [J, 0. 

We note at this point, that the terms "band gap" and 
"stop band" are related to an extended 3D PC rather 
than to finite PC samples (composed of a finite number 
of layers). However, as the thickness of a PC increases, 
band gaps and stop bands start manifesting themselves 
as high reflectivity regions. It has been shown that in 
the case of a metallic log-pile, as much as only four layers 
are already sufficient for manifestation of the band struc- 
ture features. Because of this correspondence we will con- 
tinue using the band structure terminology throughout 
the current section when discussing the high reflectivity 
(low transmittance) regions in the spectra of a finite PC. 

Let us illustrate the properties of the tungsten log-pile 
PC operating in the IR range on Fig. [JJ by FDTD sim- 
ulation. We use the same geometrical parameters, as 
in ||: the lattice constant a = 4.2/im, the rod width 
w = 1.2/xm, the rod height h = 2.4/im and the number 
of layers N = 4. The dielectric function of tungsten at 
T = 298K is used. The calculated reflectance and ab- 
sorbance spectra are shown in Fig. [TJ This PC has two 
high reflectivity regions, corresponding to the IR band 
gap and the stop band [5j, where absorption and, there- 
fore, thermal emission is suppressed, and exhibits a sharp 
absorption (emissivity) peak between the IR band gap 
and the stop band due to the interplay between the low 
light group velocity and high density of states 0, H, H, Q . 
The high emission selectivity in the near IR is a direct 
consequence of this peak. Basing on the above results, it 
was suggested 0] that by a proper scaling of the log-pile 
geometry, all spectral features could be shifted to the vis- 



ible range, thus opening the possibility of using log-piles 
as efficient light sources. 

However, our findings indicate, that high intrinsic 
losses of tungsten in the visible range wash out the rele- 
vant spectral features of the scaled log-pile geometry. We 
compare the simulation results for 4 layer perfect metal 
(no losses) and tungsten log-piles with scaled geometri- 
cal parameters (a = 0.5fim, w = 0.2/irn, h = 0.2/im) 
chosen so as to move the IR band gap edge to the vis- 
ible range. While the behaviour of the reflectance with 
the wavelength for the perfect metal log-pile is very sim- 
ilar to the one seen in Fig. [TJ the case of the tungsten 
log-pile is different. Its absorption spectrum has a single 
broad peak at the A ~ 0.5/im and falls gradually at the 
longer wavelengths, remaining quite high in the near IR 
(Fig. [5]). Correspondingly, the reflectance of the log-pile 
increases in the near IR, indicating the IR band gap for- 
mation. The stop band is completely destroyed while the 
IR band gap is strongly affected by the high absorption 
of tungsten, leading to the low emission selectivity in the 
visible range. Similar results have been obtained in the 
recent experimental study of metallic log-piles |4lj |. 

Naturally, one could think of reducing the tungsten fill- 
ing fraction in order to compensate for tungsten losses, 
which could be achieved by, e.g., decreasing the width of 
the rods. We consider the tungsten log-pile with a vary- 
ing rod width w, keeping both a and h constant. Simula- 
tions reveal, however, that, contrary to the intuition, the 
absorption increases in the near IR with decreasing the 
log width w (Fig. [3]). This behaviour can be better un- 
derstood by considering the transmission spectra of the 
perfect metal log-piles. In this case, decreasing w with 
both a and h kept constant, leads to the red shift and 
weakening of the band gaps (Fig. [3} . This can be quali- 
tatively associated with the grating-like geometry of each 
layer of the log-pile. Such a grating is characterized by 
a waveguide cutoff wavelength, corresponding to twice 
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Figure 2: Normal incidence reflection (1) and absorption (2) 
spectra of a tungsten log-pile at T = 2AQQK with JV = 4, 
a = 0.5/im, w = 0.2/im, and ft = 0.2/j.m. The reflection for 
the perfect metal PC is shown for reference (3). 
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Figure 3: Normal incidence absorption spectra of a tungsten 
log-pile at T — 2400K with N — 4, a — 0.5fim, h = 0.2/xm, 
and varying w. 

the air opening of the grating A c = 2{a — w) 8]. The 
transmittance is strongly attenuated for A > A c , indicat- 
ing the band gap formation. Decreasing w leads to the 
increasing of the air opening along with the filling fac- 
tor decrease. This results in both the red-shift and the 
weakening of the gap, allowing more radiation inside the 
structure due to the higher transmittance. This, in turn, 
leads to the higher absorption in the near IR range in 
case of the lossy tungsten log-pile, which is exactly what 
we have observed in (Fig. [3j . 

In summary, the intrinsic tungsten absorption strongly 
affects the spectral characteristics of the log-piles, ham- 
pering their application as light sources. 

4. EMISSIVE PROPERTIES OF OPAL PC 

We now turn to discussing the emissive properties of 
opal-like PCs. We consider a PC slab consisting of N lay- 
ers of spheres arranged in the FCC lattice, stacked along 
(111) crystallographic direction. Each layer of this struc- 
ture is a 2D triangular lattice with the period equal to 
a (FCC lattice constant is \/2a). There are two possible 
opal geometries: direct (tungsten spheres in air) and in- 
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Figure 4: Normal incidence transmission spectra of a perfect 
metal log-pile with N = 4, a = 0.5/im, h = 0.2/im, and 
varying w. 



verse (air spheres in tungsten slab of the finite thickness) 
opals. Here air should be replaced in practice by a refrac- 
tory supporting of filling host material, but we start from 
numerically studying the effect of a tungsten PC itself. A 
low tungsten filling fraction can be achieved with the di- 
rect opal geometry, since the sphere radius can be made 
as small as necessary. As we show in the current sec- 
tion, low filling fraction leads to higher efficiency in the 
visible. As an illustration, we present calculated absorp- 
tion spectra for direct FCC opals with filling fractions of 
2.2% and 25% and inverse FCC opal with filling fraction 
of 25% (Fig. [5]) . The absorption spectrum of the direct 
opal with / — 2.2% is characterized by a pronounced 
peak at 500nm, with amplitude several times greater than 
the absorption in the IR. The high absorbance ranges are 
also present in the spectra of both direct and inverse opals 
with / = 25%, but they are broadened with respect to the 
low filling fraction case, and their maximum amplitude 
exceeds the IR absorption by only a factor of 1.5. Note 
that in inverse opals with / < 25% the necks between the 
air cavities partially vanish. At high temperatures they 
are likely to sinter which results in thermal instability of 
this structure. This makes inverse opals less favorable for 
light applications. In view of the above arguments, we 
confine with studying the emissivity properties of direct 
opals with the low filling fraction. Below we study how 
their emissivity depend on the lattice period a, radius R 
and number of layers N. 

In Fig. [5] we present the calculated absorption for a 
two-layered direct opal PC. while scaling its geometrical 
parameters. There is a pronounced peak in the absorp- 
tion spectra, scaled linearly with a (if R/a is constant). 
By adjusting the a value, one can tune the absorption 
peak position into the visible range. The nature of this 
peak is related to the resonant absorption of the electro- 
magnetic waves in the 2D triangular layers constituting 
the opal PC. At the normal incidence, the resonant ab- 
sorption condition is fulfilled for a wavelength A when 
2ir/ A is close to the magnitude of any of the translation 
vectors in the 2D reciprocal lattice. For this wavelength 
one of the allowed non-evanescent scattered wave direc- 
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Figure 5: The absorption spectra of direct and inverse FCC 
opals at normal incidence and varying filling fraction /. The 
lattice period a = 0.45 fim, number of layers N = 2, T = 
2400K. The host medium is vacuum. 
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Figure 7: The absorption spectra of PC at different angles of 
incidence of an s-polarized wave, a — 0.45 /im, R — 0.07 \im 
and N — 2, T — 2400 A". The plane of incidence is determined 
by a normal to the layer and a 2D lattice primitive vector. 
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Figure 6: The absorption spectra of PC at normal incidence 
for varying lattice constant a, and R = 0.156a. N — 2, T — 
2400K. 



Figure 8: The absorption spectra of PC at normal incidence 
and varying sphere radius R. a = 0.45 fim, N — 2, T = 
2400A. 



tions is almost parallel to the plane of 2D layers. The 
peak, observed in Fig. HI is related to the smallest re- 
ciprocal translation vector with the length of 47ra/\/3. 
This value corresponds to A = y/3a/2, which determines 
the left edge of the peak. Peak amplitude is almost the 
same for different a values due to the weak frequency de- 
pendence of the imaginary part of the tungsten dielectric 
function in the visible. 

Note that the critical wavelength value, corresponding 
to the diffraction edge, increases with the incident angle 
value. As a result, the peak shifts towards the IR range 
while the angle of incidence increases (Fig. [7J|. 

We calculate absorption spectra for different sphere 
radii R and the number of layers N (Figs. [H [9j. One 
can see that peak amplitude grows with the increase of 
R and N, due to the larger amount of tungsten in the 
PC. It leads to the higher energy emitted in the visible 
J V is but lower efficiency E — J v - ls /Jtot- 

We now turn to the efficiency calculations for different 
values of the lattice constant 0.3/im < a < 0.7/im, sphere 
radii 0.02/tm < R < 0.2/im, and the number of layers 1 < 
N < 4. The radiation flux is obtained by integrating the 
absorption spectra over angles of incidence in the range 
0° < 6 < 90°, 0° < <j)< 360° according to Q. According 
to our calculations, for geometry with a — 0.45 /im, R = 



0.07 /im, and N = 2 (filling fraction / rj 2.2%) efficiency 
reaches the maximal value (E 15%) provided that J v j s 
is not less than J v - ls ^/2 (J v - ls ~ 15kWt/m 2 ). As has been 
mentioned above, for the bare tungsten substrate Ew = 
5%, and J v i s ,w = 30kWt/m 2 . Therefore, the efficiency 
can be increased three times with J v i s lowered by a factor 
of 2 (due to using lower tungsten filling fraction) , by using 
PC instead of the bulk tungsten. 




Figure 9: The absorption spectra of PC at normal incidence 
and varying number of layers N. a — 0.45 /im, R = 0.07 /tm, 
T = 2400K. 
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Figure 10: The top view of the PC sample with the triangular 
lattice period a — 0.55 \im. 




Figure 11: The side view of the PC sample with the triangular 
lattice period a — 0.55 fim. 



5. PC DISORDER 

In order to support the above results experimentally, a 
PC monolayer sample was fabricated on top of the quartz 
substrate. The sample consists of tungsten spherical seg- 
ment caps on top of cylindrical quartz posts arranged in 
the triangular lattice with a — 0.55 /im. The cap base ra- 
dius is 0.1 fim, the cap height is 0.02 fim, the quartz posts 
radius is 0.17 \im and the height is 0.15 fim (Figs.fTO lfTT]) . 
We choose monolayer structure since it is relatively sim- 
ple to fabricate, and the peak in the absorption spectra 
is present even if JV = 1 (Fig. [9|). The experiment has 
been carried out at the room temperature. 

In Fig. [12] we plot the comparison between the mea- 
sured absorption of the fabricated PC sample and the 
calculated absorption of a) the PC monolayer consist- 
ing of tungsten spheres with R = 0.1 [im, and b) the 
PC monolayer, consisting of tungsten caps on top of the 
quartz posts with dimensions corresponding to the ex- 
perimental sample. One can see that the peak, predicted 
by numerical calculations, is observed experimentally. In 
the case of the caps-on-posts PC the calculated absorp- 
tion peak is smaller than in the case of spheres PC, due to 
the lower tungsten filling fraction in the former case. The 
even lower peak observed experimentally is caused by the 
lattice disorder and size dispersion of the PC elements. 

One can see from Fig[T0] that experimental sample is 
characterized by some level of deviation of position and 
size of PC elements from their "perfect" values. In or- 
der to reproduce the results for such type of " imperfect" 
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Figure 12: The tungsten PC monolayer absorption spectrum 
at normal incidence at the room temperature (T = 298A"). 
The PC parameters: a — 0.55 /jm, sphere segments radius 
0.1 yitm, and height 0.02 /im. The plotted curves correspond 
to: measurement; calculation for the perfect opal monolayer 
with R — 0.1 /im; calculation for the PC with the caps-on- 
posts geometry of elements. The perfectly periodic lattice is 
assumed in calculations. 
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Figure 13: The tungsten PC monolayer absorption spectrum 
at normal incidence at the room temperature (T = 298K). 
The PC parameters: a = 0.55 /jm, sphere segments radius 
0.1 [im, and height 0.02 /im. Comparison of the measured 
spectrum to the calculation with both the size dispersion and 
the lattice disorder taken into account: results for the varying 
lattice disorder parameter p max are presented (see text). 



PC numerically, we use the following scheme: a suffi- 
ciently large PC area is considered, with the XY dimen- 
sions much greater than the lattice constant. Each PC 
element is shifted randomly from its lattice node, the 
displacement being characterized by a vector with abso- 
lute value p = p max ?7 p in the direction, specified by an 
angle (f> = 2~Kr\$ with respect to the X axis, where p m ax 
is the maximal possible displacement. The volume of 
each PC element can vary around its mean value Vq as 
V = Vq + AV„ lax (2r]v — 1), where AVnmx is the maxi- 
mum possible variation, and r\ P ^y are random numbers 
in [0, 1] interval. According to this scheme, an imperfect 
PC is generated by using independent random values r\ 
for assigning position and volume of each PC element ac- 
cording to the expressions given above. Generated struc- 
ture is placed into the FDTD computational cell with 
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periodic boundary conditions in the X and Y directions. 
This structure is randomly disordered at the small scale 
(~ a) , and uniform at the large scale (>■ a) that leads to 
the self-averaging of its optical properties (42J. Absorp- 
tion spectrum of this structure with appropriately cho- 
sen parameters p m ax and AVmax appears to be a good 
approximation for the absorption spectrum of the imper- 
fect experimental PC. 

In Fig. [13] we present comparison of the calculated 
absorption spectrum for generated structure with a = 
0.55/xm and experimental data. Computational cell lat- 
eral size is 8x8 PC unit cells with the period a (further 
increase of lateral size does not change the results sig- 
nificantly). The base radius of the sphere segments is 
varied, while its height is kept constant. The maximum 
volume variation is set to AVmax = 0.5Vo- The maximum 
shift p max = 0.3a gives the best agreement of the calcu- 
lated peak height with the experimental data. Note that 
the peak height decreases with the increase of parameter 
Pmax corresponding to PC lattice disorder (Fig. Q2]). At 
the same time, the IR part of the spectrum depends only 
on the size of the elements (e.g., see Fig. |8]), since the 
wavelengths in IR are much greater than the lattice con- 
stant, and PC can be described as the effective medium 
in this spectral range. 

In summary, the " imperfectness" of the experimental 
PC leads to decrease of the emission in the visible, how- 
ever, the IR emission is still small. 



6. PC IN A HOST MATRIX 

The scattering elements of a PC should be embedded 
in a refractory material (host). There is only a small 
number of refractory materials, resistant to high tem- 
peratures and applicable for this purpose (i. e., yttrium- 
stabilized zirconia and hafnia). However, when being 
heated, these materials emit in the IR. We find that it 
leads to a substantial efficiency reduction of the whole 
structure (PC and host). In this section we illustrate 
this result using an example of a tungsten two-layer PC 
embedded in the hafnia (Hf02) host (Fig. [14J inset) with 
the width d: 



d = 2R+{l + l/VS)a. 



(3) 



We start with a non-absorbing host, taking the hafnia 
dielectric permittivity at room temperature (e = 4.41). 
To avoid oscillations of the absorption curve caused by 
the finite thickness d of the host slab, we consider the PC 
embedded in an infinite environment with s = 4.41. The 
change of the environment dielectric function e from 1 to 
4.41 leads to a red-shift of the absorption spectrum by 
a factor of y/e. It is explained by the ^/e-fold reduction 
of the wavelength in the external environment (compare 
blue dotted curves in Figj7] and Fig. [T4] corresponding 
to the same PC geometry in vacuum and non- absorbing 
hafnia environment). As a result, optimal parameters for 
the PC in host can be obtained by scaling down by y/e 




0.6 0.8 1 

wavelength (urn) 

Figure 14: The absorption spectra of the two-layer PC at 
normal incidence (T = 24007^) . 1 - PC is embedded in hafnia 
at the room temperature, a — 0.45 /im, R — 0.07 fim; 2 - 
the same as in 1, with a = 0.214 \im, R = 0.033 )im; 3 - the 
finite thickness of the host is taken into account, d is given 
by @. Inset: Two- layer PC embedded in plane-parallel host 
with width d. 



that ones for vacuum environment (see green curve in 
Fig. ED. 

The finite thickness d of the embedding host leads to 
the interference of the waves multiply reflected from the 
host- vacuum interfaces. It leads to additional Fabry- 
Perot oscillations in the absorption spectrum. However, 
it does not affect the location of the main absorption 
peak. It can be seen in Fig. [TH where the results for the 
PC in the infinite external environment (green curve) and 
in the finite slab of width d defined by © (red dotted 
curve) are compared. 

Hafnium becomes strongly absorbing at high temper- 
atures. Its dielectric function at T = is well de- 
scribed by Lorentz approximation: 



e(u>) 



1 



Aeuj% 



ujq — 2iujj — to 



2 ' 



(4) 



Here Ae = 3, w = 1.2 • 10 16 rad/sec, 7 = 1.0 • 
10 16 rad/sec. 

Effect of the host absorption can be examined by as- 
suming 7 = s ■ 10 16 rad/sec with parameter s increasing 
gradually from (zero absorption) to 1 (value from ex- 
pression |4} . One can see from Fig. [T5l that for the small 
s values the absorption spectrum is mainly determined 
by the tungsten spheres. Increasing s leads to the higher 
hafnia absorption in the IR and a significant decrease of 
the efficiency of the whole structure (7% for s = 1). In 
summary, the photonic absorption peak is washed out 
by the absorption in the host matrix, necessary for em- 
bedding PC elements, so that a practical realization of 
PC-based light sources becomes difficult. 



7. CONCLUSIONS 

In this paper, we studied the possibility of using PCs 
as high-efficiency light sources. We considered different 
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1 r 




0' 

0.4 0.6 0.8 1 1.2 

wavelength (liiti) 



Figure 15: The absorption spectra of the PC at normal in- 
cidence (T = 2400K). The PC parameters: a = 0.214 fim, 
R = 0.033 /J,m, N = 2, d is given by ([3]). Different values of 
the damping coefficient in the Lorentz formula for absorbing 
hafnia are considered (see text). The same FDTD parameters 
as on Fig. [6] are used. 

PC geometries (log-piles, direct and inverse opals) and 
showed that a direct opal geometry is more appropriate 



for our purpose. By changing geometrical parameters of 
a direct opal PC, one can achieve high selectivity of emis- 
sion in the visible range. In real PC samples, a lattice 
disorder and a size dispersion of the constituent elements 
are present. This leads to a decrease of the peak emissiv- 
ity in the visible, however, emissivity in the IR remains 
small. 

For a practical application, the direct opal lattice needs 
to be embedded in a refractory host matrix. Known re- 
fractory materials are characterized by a high emissivity 
at the temperatures of interest. This leads to a lower 
PC peak emissivity in the visible. Therefore, implement- 
ing PC based light sources requires that low- absorbing 
refractory materials should be explored. 
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